Short time growth of a KPZ interface with fiat initial conditions 
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The short time behavior of the 1 + 1 dimensional KPZ growth equation with a flat initial condition 
is obtained from the exact expressions of the moments of the partition function of a directed polymer 
with one endpoint free and the other fixed. From these expressions, the short time expansions of 
the lowest cumulants of the KPZ height field are exactly derived. The results for these two classes 
of cumulants are checked in high precision lattice numerical simulations. The short time limit 
considered here is relevant for the study of the interface growth in the large diffusivity/weak noise 
limit, and describes the universal crossover between the Edwards- Wilkinson and KPZ universality 
classes for an initially flat interface. 

PACS numbers: 05.40.-a, 05.20.-y, 05.70.Np 



I. INTRODUCTION 

The growth of interfaces in the presence of noise can 
be classified into several universality classes. When the 
growth rate does not depend on the slope of the inter- 
face, the growth process falls in the simplest Edwards- 
Wilkinson (EW) class. When instead the growth rate is 
slope dependent (e.g. in the so-called lateral growth), 
the process falls into the class defined by the non-linear 
continuum Kardar-Parisi-Zhang (KPZ) equation [TJ [5] 

d t h = vV 2 h+±\ (Vh) 2 +£(x,t), (1) 

where h(x,t) is the interface height, v the diffusivity, Ao 
the strength of the slope dependent growth (with Ao = 
giving the EW model). £(x,t) is the stochastic noise, 
chosen centered gaussian with short range correlations 

£(x, t)£{x',t') = R 6 (x- x')5(t - t') . (2) 

Concerning the space dependence, the most usual choice 
is to take uncorrelated random disorder, i.e. R^(x — x') — 
DS(x-x'). 

In one spatial dimension, the KPZ universality class 
shows a fairly robust anomalous scaling exponent [3] for 
the fluctuations of the height of the interface h(x, t) ~ 
t 1 / 3 and indeed such anomalous behavior at large time 
has been proved for several discrete solvable models [3HZ] 
which are believed to belong to the KPZ class. However, 
this exponent is only one of the facets of the universal- 
ity of the KPZ equation: further universal information 
is encoded in the full probability distribution function 
(PDF) of these fluctuations, but their exact calculation 
is an extremely difficult task which is complicated by the 
fact that, even after a long time, the system keeps some 
memory of the initial conditions. Remarkably, these ini- 
tial conditions can be classified in a few subclasses, each 
leading to a distinct universal result for the statistics of 
the height field at large time [9] . An impressive the- 
oretical progress has been recently achieved and has led 



to exact solutions directly on the continuum KPZ equa- 
tion for the wedge (or droplet) [TUHT3] . flat E>] and 
stationary p~6] initial conditions. In the first two cases 
the PDF of the height h(x,t) at a given point converges 
at large time to the so-called Tracy- Widom GUE (Gaus- 
sian unitary ensemble) and GOE (Gaussian orthogonal 
ensemble) universal distributions [17j , for droplet and flat 
initial conditions respectively. Further impetus to the 
field has been given by recent experiments on turbulent 
liquid crystals [HI [19] in which these two long-time pre- 
dictions have been confirmed with high accuracy. 

In the literature, much emphasis has been given to the 
long time limit, mainly because of the connection with 
random matrix theory. However, the above mentioned 
exact solutions for the KPZ height distributions are valid 
for arbitrary times. Indeed these solutions can be ex- 
pressed in terms of Fredholm determinants with rather 
complicated kernels, from which it is not always easy to 
extract the limiting behavior for long and short times. 
It is then interesting to obtain, by simpler means, the 
small time behavior in an explicit form, and to confirm 
it in numerical simulations. This has been achieved in 
the case of the droplet initial conditions [13 [5D] , and the 
aim of this paper is to present a similar result in the case 
of the flat initial condition. As discussed in more details 
below, there are generically three time regimes: 

(i) a non-universal very short time regime t ~ tt 
where the growth depends on short scale details 
(e.g. small deviations from flat initial condition, 
precise form of R^(x), etc. . . ); 

(ii) a short time regime tf -C t <C t* where the 
crossover from EW to KPZ takes place; 

(iii) a large time regime t ^> t* where KPZ scaling 
holds. 

The short time regime (ii) is of particular interest in the 
high diffusivity/weak noise limit when tf <C t*\ it occurs 
in a broad window of time scales and it is independent of 



2 



short scale details. It is this universal crossover that we 
are studying here. 

The paper is organized as follows. In the next section 
we discuss the mapping of the KPZ equation to the di- 
rected polymer. In Sec. |III[ we report the small time 
expansion of the moments of the partition function of 
the DP obtained in Ref. [15] and we check them in nu- 
merical simulations. In Sec. |IV| we calculate analytically 
the small time expansion of the connected moments of 
the height held and in Sec. [V]we check them by numer- 
ical simulations. Three appendices contain some more 
technical calculations. 



II. MAPPING TO THE DIRECTED POLYMER 



Via the Cole-Hopf transformation, the KPZ equation 
can be mapped onto the directed polymer in a random 
environment which is an equilibrium statistical physics 
problem [I] [5T1 [25]. A growth starting from a droplet 
initial condition is mapped onto a hxed endpoints poly- 
mer, while a flat initial surface translates to a directed 
polymer with one endpoint fixed and the other free [15j . 
Indeed the canonical partition function of a directed poly- 
mer x(t) at temperature T in a random environment is 
defined in the continuum by the path integral 

Z(x,t\y,0) = [ X{t) ~ X Dxe-± fo^lhm 2 +v( x (r),T)] ^ (3) 

Jx{0)=y 

and maps to the KPZ equation after the identifications 
A 



2v 



h = UZ, 2v = T, M£{x,t) = -V(x,t). (4) 



A Gaussian noise £(x, t) corresponds to a random poten- 
tial V(x,t) which is a centered Gaussian with correla- 
tor V(x,t)V(x',t) = R v (x - x')5{t - t') with R v (x) = 
XqR^(x). The white noise in KPZ corresponds in poly- 
mer language to disorder with ^-correlations 



V(x, t)V(x', t) = cS(t - t')8{x - x'), 



= D\%. (5) 



This mapping is valid in the bulk and does not depend on 
the KPZ initial condition which translates into conditions 
for the endpoints of the polymer. For the KPZ equation 
with flat initial condition, one should consider the parti- 
tion sum with one fixed endpoint (at x) and another free 
(at y) [TJJ [TS] resulting in the partition function 



Z(x,t) 



dyZ(x,t\y,0) . 



(6) 



The recent analytical progress has been made possi- 
ble by the calculation of the moments Z(x,t) n of the 
DP partition sum. By replicating the partition function 
Z(x,t), the DP is mapped (23] onto the quantum me- 
chanics of a bosonic system of n particles interacting with 
an attractive delta-function potential, i.e. the celebrated 



Lieb-Liniger model [53] . This model is integrable via the 
Bethe Ansatz and the eigenstates are known both for re- 
pulsive [53] and attractive interactions [55] which is the 
case of our interest. The moments can be expressed as 
a sum over these eigenstates [TU [T5] (generically labeled 
by [i in the following) 



z{x,ty 
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in terms of the many-body wave-function 4 1 ^(j/i, • • • , y n ) 
and of the eigenenergies of the state /x. In the infinite 
system the eigenstates are easily enumerated, being or- 
ganized in clusters of bound particles, called strings. The 
norms of the states ||/i|| and the equal points wave func- 
tions have simple expressions [26] and led to the time- 
dependent PDF starting from a droplet initial condition 
[TTT 115] . The integral over the £/, in Eq. Q , necessary 
to treat the flat initial condition, is more delicate but 
was handled in Refs. [HI [TS] leading to the moments 
Z{x,t) n for arbitrary n. From these the moment's gen- 
crating function at all times has been written in terms of 
a Fredholm Pfaffian [T31 ITS] (the square root of a Fred- 
holm determinant). This allowed to prove that the PDF 
of hiZ(x,t), i.e. of the height field h(x,t), converges at 
large times to the GOE Tracy- Widom distribution. 

Here we follow a different route. We recall in the next 
section the exact expressions for the lowest moments n = 
2, 3, 4 and from them we extract the small time cumulants 
of In Z, i.e. of the KPZ height field. 



III. MOMENTS Z n AND THEIR SMALL TIME 
BEHAVIOR 

For flat initial condition, the one point distribution of 
Z(x, t) does not depend on x because of translational 
invariance. Thus in the following, we simply denote 



Z = Z(x,t). 



(8) 



Since we are dealing with an initially flat interface we 
must have Z n = 1 at t — (which is a non-trivial condi- 
tion in terms of the Bethe Ansatz). Taking the average 
of Eq. j3j over the Gaussian disorder gives the mean 
partition function 



Z(x,t\y,0) = 



1 (x-y) 2 MyjO) . 

e ~ art , (9) 



V2irTt 

and so from the integral in Eq. ^ we have 



s" ', with 



Rv(0) 
2T 2 ' 



(10) 



To eliminate this non-universal (self-energy) contribu- 
tion, it is convenient to define 



Z/Z, 



(11) 
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that by construction satisfies z = 1 at all times [28 . This 
will be useful later for comparison with the numerical 
simulation of lattice models. 

All results for the continuum DP and KPZ models are 
expressed in terms of a dimensionless parameter 



1/3 



with t* = 



2T 5 2{2vf 



d 2 K 



(12) 



where, in the language of the DP, t* is the crossover 
time scale between the Brownian diffusion at small time 
(i.e. A < 1) and the glassy large time behavior (i.e. 
A > 1). Within the context of the growth model, t* 
is the crossover scale between Edwards- Wilkinson and 
KPZ regimes. Note that a spatial crossover scale can be 
also defined as x* = ^fut* = T 3 /c = (2^) 3 /£> 2 Ag. Both 
scales become very large in the large diffusivity limit or, 
equivalently, in the weak noise limit. 

It is important to recall that for any given microscop- 
ical model with a cutoff (e.g. a lattice model) there are 
additional time/space scales. The easiest example is the 
same continuum KPZ equation (or DP model) with a 
disorder correlated over a non-zero correlation length r*y, 
i.e. R^(x) is a function decaying on a scale rf. Then, 
it is easily shown [TT] that if x* ^> r/ one can replace 
R^{x) —> DS(x) in which D = J dxR^(x). More gen- 
erally, the condition for the existence of the universal 
short time regime studied here is that t* and x* must 
be much larger than any characteristic microscopic scale 
-generically called r/ and tf here- such as the lattice 
spacing for a lattice model. Note also that if the initial 
condition is not perfectly flat on scales of the order of Tt, 
this will also not affect any result as long as x* >• r/ . Of 
course, the very short time/space regime with t <tf and 
x <rf is non- universal. 

We now recall the results of Ref. [15] for the four lowest 
moments, together with their small time (i.e. small A) 
behavior 



z 2 = e 2A [l + erf(A 3 / 2 v / 2)] 



(13) 



= l + 2\/ y A 3 / 2 + 2A 3 + ^/ T A 9 / 2 + 0(A 6 ), 

V 7T 3 V 7T 



with [29 erf (z) = 4» /* dte 



4e 8A3 - 2e 2A3 - 2e 8A3 erfc(A 3 / 2 2\/2 
+e 2A3 erfc(A 3 / 2 y2) 



(14) 



1 + 6 J ^A 3 / 2 + 14A 3 + 40W ^A 9 / 2 + 0(A n / 2 ) 




0.001 



FIG. 1: From top to bottom the moments z 3 — 1 (solid line, 
squares) and z 2 — 1 (solid line, circles) for different values 
of A. So lid lines c orre spond to the analytical predictions in 
Eq. (14 1 and Eq. (131. Averages are performed over 15 ■ 10 6 
512 and c 



samples of size t 
parameters. 



1. There are no adjustable 



with erfc(x) = 1 — erf(x), and finally 



„8A a 



4e 20A3 erfc(3V2A 3 / 2 ) (15) 



-4e 8A3 [-2erfc(2A 3 / 2 ) + e 12A3 erfc(4A 3 / 2 )] 



/ 2x+l 



x+1 ' C 



N /, dX 4r.-(4.r(,r + .S) + .",) 

7T 



l + 12^A 3 / 2 + (44+^)A 3 
+8(21 + 8^2) ^|a 9 / 2 + 0(A 6 ) , (16) 



where the series expansion of the integral is performed in 
Appendix [5] 

Before embarking in the calculation of the cumulants of 
In Z, we now report the results of numerical simulations 
for the determination of z n for n = 2, 3. As explained 
in more details in Sec. [VJ the simulations are performed 
for a directed polymer on a square lattice. We also con- 
sider the high temperature limit which ensures that all 
details of the lattice become irrelevant and the results 
can be expressed as functions of the single parameter A. 
The procedure and the identification of A on the lattice 
have been introduced already in [TTJ [27] and is described 
again in Sec. [V] We report the numerical data for z n in 
Fig. [T] which are found to be in excellent agreement with 
our analytic predictions up to A sa 0.6, while some devia- 
tions at larger A are evident. These deviations are caused 
by the undersampling due to the growing importance of 
the tails in the distribution of z, and will be properly 
explained in Sec. [V] 
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IV. CUMULANTS OF In Z AT SMALL TIME 

From the above formulas for z n and following the pro- 
cedure described in Appendix [Bj we obtain the small A 
(i.e. small time) expansion of the first four cumulants of 
the free energy 



In z 



(In; 



(lnz) 3 



(lnz) 



A 3/2 



6 



A 



3 7T 

U»v / 2--)a/ ¥ A 9 / 2 +0(A 6 ), 

7T / V 7T 



(17) 



^2 - 



= 2 



7T 

/106 - 40 

AT 

/|A3/ 2+ (^- 6 )A3 (18) 

(19) 



(248 - 96^2 - ^) ^|a 9 / 2 + 0(A 6 ), 



64^2 + — -192 



A 9/2 + 0(A 6 ), (20) 




FIG. 2: Sketch of the directed polymer model analyzed in 
numerical simulations. The blue solid line corresponds to a 
polymer growing over the square lattice under the hard con- 
strain condition. 



and of coursefln Z) p = (In z)p for p > 2. As explained 
in Appendix [bI in order to compute the next term 0(A 6 ) 
in the small time expansion, or the fifth and higher cu- 
mulants, we would need the fifth moment z 5 that we 



did not analyze here, but which is in principle known 
[15] . A simple check can be performed on these formu- 
las, namely one can compute the series expansion in A 3 / 2 
of exp(^ p=1 ^(lnz) p ) and check that the expansion of 
all z n for n — 1,2,3,4 given above is reproduced up to 
order o(A 9 / 2 ). Although this procedure also allows to de- 
rive Eqs. ( [l7p0 1 by adjusting the coefficients of the series 
expansion in A J ' 2 , the method described in Appendix |b| 
is more systematic. 

For short times, the dominant term in the PDF is the 
variance (lnz) 2 which increases as t 1 / 2 . Using InZ = 
Aq/i/(2z/), one finds 



h 2 =D 



t 

2tw 



0{\ 2 t) 



(21) 



Hence the first term of the expansion of h 2 is indepen- 
dent of Ao, the coefficient of the non-linear growth in the 
KPZ equation. It corresponds to the Edwards- Wilkinson 
Gaussian scaling regime Sh ~ i 1 / 4 , also found in Ap- 
pendix [C] cf . Eq. ( C3 1 , where we derive the leading 



short time behavior for the average height and variance 
using perturbation theory directly on the KPZ equation. 

The third and fourth cumulants behave as t and t 3 / 2 
respectively, suggesting that the fourth cumulant is sub- 
dominant and that the first corrections to the EW gaus- 
sian scaling are given by the third cumulant as 5h ~ i 1 / 3 , 
which is the form of the KPZ scaling. Indeed, it is inter- 
esting that the third cumulant is linear in t both at short 



and large times (but with an the amplitude going from 
= 0.360563 to ^ OE = 0.598268, see below). 
From the above we can compute the skewness 



71 



(lnZ) s 



. *l< (22) 

[(InZ) 2 ] 3 / 2 \fr 

7T (67 - 48^2 + 9vr) -86 



+ 



-A 9 / 4 + 0(A 15 / 4 ) 



(2tt) 3 / 4 

= 0.178865 . . . A 3/4 + 0.0138 . . . X 9/i + 0(A 15/4 ), 
which at short times scales as 71 ~ i 1 / 4 and the kurtosis 



72 



(ln^) 4 
[(InZ) 2 ' 



(40 - 24tt + 8V2it) J - A 3/2 + <3(A 3 ) 
= 0.115565... A 3/2 + 0(A 3 ), (23) 



which scales as t 1 / 2 . 

Now we recall that at large time one can write [21 H"5] 



Xoh 
2v 



= InZ = Voot + Xr] t 



(24) 



such that rjt converges to the GOE Tracy- Widom distri- 
bution lim Prob(?7 t < s) — i*i(s). The skewness and 

t— >oo 

kurtosis thus converge for large times to their GOE val- 



71 -> 1?° E = 0.29346452408 . 

72 -> l2° E = 0.1652429384 . . 



(25) 
(26) 



consistent with a crossover for A ~ 1.6 ± 0.3. The 
amplitude of the (non-fluctuating) linear term is non 
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k k 



FIG. 3: From top to bottom, the cumulants (15- 10 6 samples) 
In 2 z (dashed line, circles), — lnz c (dashed line, triangles) and 

q C 

In z (dashed line, squares) for t = 512. The dashed lines are 



the analytical predictions in Eqs. 1 171, (181, and (19 1 up to 
0(A 9/2 ), with c= 1. There are no adjustable parameters. 



FIG. 4: Finite-size effects for In 2 z . Solid line: analytical 



prediction Eq. Q 1 8 p - Numerical data, from top to bottom 
t = 128 (green circles), t = 256 (blue circles), t = 512 (red 
circles). Averages are performed over 15 • 10 6 samples with 
c= 1. 



universal, Voo = vq — c 2 /12 = Vq — D 2 Xq/12, where 
<^v(o) _ x r^(o) j g ^ e amplitude at short time 



2T 2 



(after the very short time regime t^$> tf). Note that the 
difference — v is universal. At large time one also 
has that lnz = \nZ — \nZ = A/ii — D 2 Apt/12 is univer- 
sal, where fif OE = -1.2065335745820... is the mean 
of the TW distribution, while (\nZ) — > A 2 /i2 where 
^goe = 1.607781034581 ... is the variance of the TW 
distribution. 



V. NUMERICAL RESULTS 



To avoid numerically instabilities we divided all partition 
functions at fixed r by the biggest one and record its log- 
arithmic value. As in the model in the continuum, also 
on the lattice Z(x, t) blows up exponentially in time, as 
can be seen by averaging the sum over all possible paths 



z(t) = j2 n e ' pv{x) = 2 '^ 2t/2 



(29) 



7* xeit 



For this reason we work numerically with the ratio 
ln(z) = \n(Z/Z), which remains small, but exhibits 
strong fluctuations. 



Numerical simulations are performed for the square 
lattice model depicted in Fig. [2j Directed paths grow 
along the diagonals of the lattice with only (0, 1) or (1, 0) 
moves (hard constrain condition), starting in (0,0) and 
with the second end left free. To each site of the lattice 
is associated a i.i.d. random number V(x, t) (here we use 
a Gaussian distribution with variance equal to 1). The 
time coordinate is given by t = i + j and the space co- 
ordinate by x = (i — j)/2 (see Fig. [2|. The partition 
sum over all paths 74 growing from (0, 0) up to time t is 
defined as 



(27) 



it 



The partition function satisfies the following transfer ma- 
trix recurrence relation implemented in our simulation 



(Z x i t +Z x+ i t )e-P v *-^ 



(28) 



with Z x = 5 x q. The free end partition function is com- 
puted by summing over all endpoints Z(t) = Y] T Z(x,t). 
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FIG. 5: Finite-size effects for In 2 . Solid line: analytical 
prediction Eq. ( |17[ ). Numerical data, from top to bottom 
t = 256 (blue triangles), t = 512 (red triangles). Averages are 
performed over 15 ■ 10 6 samples with 5=1. The large error 
bars when A < 10~ 2 are due to the vanishing value of In z as 
A -> 0. 
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z 



FIG. 6: P(z) with z = | for A = 1.26,0.58,0.27,0.126 and 
0.058, from left to right. Histograms are obtained from nu- 
merical simulations with t = 512, c = 1 and 15 • 10 6 samples. 
When A is very small, P(z) is self- averaging. When A grows, 
a heavy tail is developed and z typ <C (z) = 1. 




20 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 

10 6 iterations 

FIG. 7: Estimator of ~z? - 1 given by M N = i 2? - 1, 

where z$ is the rescaled partition function of a single disorder 
realization and N is the number of realizations (iterations), 
for A = 1.26. The sudden variations around 10 • 10 6 samples 
shows that one single event contributes to a finite fraction 
(around 10%) of the whole sum. 



In the limit of high T, the statistical fluctuations of z 
only depend on the unique dimensionless variable 



_ /c 2 ^\i/3 



(30) 



which is the lattice version of Eq. (12 1. In the high 



temperature regime, the parameters k and c can be com- 
puted explicitly [27]. Indeed c 2 is just the variance of 
the uncorrelated random numbers. Instead k can be 
extracted from the model without disorder, for which 
the polymer behaves like a particle diffusing on a one- 
dimensional lattice (x being the particle position at time 
t). The mean square displacement of the particle is given 
by (x 2 (t))x — tx/n. Within the normalization used in 
this paper, we have k — AT. 

Using this algorithm, we have numerically computed 
the cumulants on the square lattice at high temperature 
and compared them with the analytic predictions in Eqs. 
([17]), and The data for In" / for n = 1, 2, 3 are 
reported in Fig. |3J The agreement with the analytical 
predictions is excellent, which is even more impressive 
when we consider that these figures are produced without 
any fit parameter. In Figs. [4] and [5] we show in more 
details the convergence to the theoretical value for a fixed 
value of A as a function of polymer length. The increase 
of the polymer length t is equivalent to heating up the 
system, hence approaching the universal prediction of the 
high temperature regime. 

The analytic predictions for the moments of z are exact 
for all A. However, we can see in Fig. [T] that precision 
is quickly lost above the threshold A ~ 1. This is due 
to the fact that, for large A, typical values of z strongly 
differ from the average value z — 1. The moments of z 
are then dominated by rare occurrences of very large z 
induced by the presence of heavy tails. This is shown 



in Fig. [6] where we see that, as A grows, the mode of 
the distribution quickly goes to while the tail becomes 
fatter. A simple example of this peculiar behavior lays 
in the log-normal probability distribution which is char- 
acterized by an exponentially small typical value and a 
i 2 

heavy tail ~ e n z . Here, for large A, the heavy tail be- 



haves as e" 



-In 3 / 2 ; 



with the exponent 3/2 corresponding 



to the Tracy- Widom asymptotic behavior. In practice, 
because of this tail, the moments estimators converge 
extremely slowly, even for important sampling. An ex- 
ample is shown in Fig. [7J where we estimate z 2 — 1 for 
A = 1.26 with N = 15_- 10 6 samples. While the predic- 
tion from Eq. ( 18 ) is z 2 — 1 = 109.023, we found a value 



around 60 from numerical simulations. This discrepancy 
is explained by the central limit theorem that predicts 

we see 

14 



20 



fluctuations of order ~ (z 4 /N) 1 / 2 . Using Eq 
that z A grows very fast in A and would require iV = 10 
samples to have a good estimation of z 2 for A = 1.26. 



VI. CONCLUSION 

In this paper we have studied the KPZ equation 
with flat initial conditions and extracted from the re- 
sults of Ref. [H] the short time behavior of the con- 
nected moments of the distribution of the height field 
at a given point. The importance of these results stems 
from the proof of the existence of a short time universal 
regime which describes the crossover from the Edward- 
Wilkinson to the KPZ growth and which can be observed 
when the diffusivity is large or the noise is weak. We have 
compared our analytical predictions to high precision nu- 
merical simulations of a discrete model, which shows how 
this universality arises. A part from the per se theoretical 
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interest, these predictions should be useful also in future 
experiments in which the parameters of the growth could 
be varied and controlled in a more refined way so to easily 
access this universal crossover. 



Appendix A: Expansion of an integral 



Thus, putting together the three pieces, we have 



7=l-12W-A 3 / 2 + (44+-)A 3 

V IT 7T 



-8(21 + 8\/2)^A 9/2 + 0(A 6 ). (A8) 



We need to compute the small A expansion of the in- 
tegral 



-8\ 3 x 



/>oo 

I = dxf(x)t 
Jo 

where f(x) and its large x expansion are 



(Al) 



f(x) = 48 



2x+l _ V2x+1 
\/4x + l x+1 

4tt(4x(x + 3) + 5) 



3 x s /2 _ 3 (21 + 8 ^ a .-6/2 + Q(a; -7/2 )- (A2) 

87T 



7T 

It is convenient to write 
I 



= / dxf(x) + / 
Jo Jo 



<l.rf (.,■)+ I dx-^(e^ x3x -1) 



+ 



/>00 

/ dxfx{x){e 
Jo 



1), (A3) 



with fi(x) = f(x) — 3x 3 / 2 /n. Two integrals are easily 
done, giving 



I = 1 



^^A 3 / 2 + dxh{x){e-* xi * - 1), (A4) 

and the remaining integral is 0(A 3 ). This can again be 
written as 

/>oo />oo 

/ dxh{x){e^* - 1) - (-8A 3 ) / dxxh{x) 
Jo Jo 

+ / dxhix)^- 8 * 3 * -l + 8\ 3 x) 
Jo 

/ 9zL\ /*°° 

= (44 + — J A 3 + / dxh{x){e- 8x3x - 1 + 8A 3 x), 

(A5) 

where the remaining integral is now 0(A 9//2 ) and can be 
split again as 



dxh{x){e- &x3x - l + 8A 3 a;) = 



dx 



3(21 



,-5/2 



8tt 

-8A 3 a; 



(e 



+ / dxf 2 (x)(e 
Jo 

with 

f2(x) = h(x) 



1 + 8A 3 x), 



1 + 8A 3 a;) 
(A6) 



3(21 + 8^/2) ... 
-i — - — —^x~^ 2 . (A7) 
on 



Appendix B: From the moments of Z to the 
moments of In Z 

In general the knowledge of the moments z n for some 
low integer n does not allow to extract much information 
about the cumulants (log z) n . In the present case, how- 
ever, for small time (small A), z is concentrated around 
its mean value z = 1 and this allows to obtain the be- 
havior of the cumulants at small times. 

Let us write z = 1 + u with u = and compute its 
connected moments. In order to lighten the notation we 
introduce the notation /i p = zP c . Using the expressions 
for z nC in the main text, u™ c are given by 

^2 = u 2 = z 2 — 1 

= 2\/^A 3 / 2 + 2A 3 + |\/^A 9 / 2 + 0(A 6 ), (Bl) 

V 7T 3 V 7T 

^ 3 = ^= -3(^- 1) 

= 8A 3 + 32^A 9 / 2 + (9(A 6 ), (B2) 
H 4 = «4 C =?-1-4(F-1) + 6(F-1)-3(F-1) 2 

= 64%/2^A 9 / 2 + 0(A 6 ). (B3) 

Given the above trend it is reasonable to assume that 
z p c = uP c = 0((A 3 / 2 ) P_1 ). Based on this assumption, 
we want to construct a systematic series expansion of 
(log 2)™ in powers of the cumulants of z. The reasoning 
is the following. First we write 



7" c 

^2 ^y(log^)" C = InF = In (1 + • 



n=l 



In 



00 



rjr- l)..(r-fc + l) — 



k=i 



k\ 



). (B4) 



Expanding the rhs in powers of r, we obtain formally each 
(log z) n as an (infinite) scries of the moments u k . The 
moments u k can themselves be expressed as functions of 
the cumulants \i p by writing 



fc=2 



to! 



exp 



p=2 F 



(B5) 



and identifying order by order in w. We can now replace 
/i p — > a v ~ x \i v where a is to be set to unity at the end. 
This replacement allows to keep track of the order in A 3 / 2 



of each cumulant. Using Mathematica, it is now simple 
to first generate the expansion ( B5 ) , truncate it to a given 



order in a, and then insert it in Eq. (B4|. During this 
process, we see that e.g. it 4 = 0(a 2 ) = u 3 , i.e. in Eq. 
one must keep a few more orders compared to Eq. 
Since we have not computed z 5 = 0(a 4 ), we can 
only get our cumulants up to order a 3 . Doing so we 
obtain 



B4 
B5 



log (z) = - 



log(z) 



log(z 



|3 - 



log(z) 4 



a^2 . 2 /A*3 3/x| 



+a 



2 V 3 

if 5 M ? 



(B6) 



afi2 + a 



f + 2^2 -^)+0(a% 

'5^ 2 

~2 



(B7) 



V 3 



12 / 



a 2 [1^3 



3 M l) (B8) 
+a 3 ( - 22^ + 15^2 - 3 |i) + 0(a 4 ), 
a 3 (20/4 - 12^3^2 + A**) + 0(a 4 ). (B9) 



simple algebra gives 



2^ 



dt\e 



(C3) 



2^* (g) 2, g 2 



," V 2tt' 

2 



where the last equation is valid for t ^> H/i/, i.e. 
away from the (non-universal) very short time regime. 
Here D = R%(q = 0) = / dxR(x) is the only mem- 
ory of the short scale details and thus for t ^> £/ one 
can set Rf(q = 0) = D, i.e. a delta-correlator in 
space. Using the correspondence InZ = Xoh/(2v) and 
A 3 / 2 = £>A^(t/8) 1/2 /(2^) 5 / 2 one recovers exactly the 
leading term in Eq. (21 1. 



The discussion of the average height h is more subtle 
because we need to retain R^(x) in an essential way, as 
there are non-universal contributions. For a flat initial 
condition, h(x,t) is x independent, hence the following 
equation is exact at all times 



An 



(C4) 



Setting a = 1 and replacing the [i v by their actual values 
above, we find the result given in the text. 



Appendix C: Short time perturbation theory for the 
KPZ equation 

As a further final check, we recover here the leading 
short time behavior for the first and second cumulants 
of the height directly from the perturbative expansion of 
the KPZ equation. We start from the second cumulant 
which is easier. The KPZ equation can be studied in 
perturbation in An which is equivalent to short time. This 
is clear from the definition of A in Eq. ( 12 ) which gives 
the perturbative parameter A 3 / 2 oc y/t/t* dc v^Aq. We 
can write h = + + ... where = O(Aft). 



With the flat initial condition, the lowest order is just 
the Edwards- Wilkinson result that in Fourier space is 



h {0) - 



(CI) 



which leads to the variance 



' l q,t' l q\t — 



ft 

2ir6(q + q') / dhe' 2 ^ 1 R $ (q), ( C2 ) 
Jo 



where R${q) is the Fourier transform of the noise corre- 
lator R(x), assumed to be of range r/ in space. Then 



At small time we can use 



A 



ftfc = ^(v*«»)3 = g J ikbW -*-*"**)■ (C5) 

Splitting this term in two pieces and integrating each of 
them separately over time, we obtain 



h= — ; 1 - — I —R^q) — — 3 . (C6) 



Av / 2tt 



One recognizes the same integral entering the second mo- 
ment, and so for t -C t* we have 



v t ■ 



4^ ' 



(C7) 



which indeed reproduces exactly, for t ^> tj, the lead- 
ing negative correction in Eq. (17). The first term lin- 



ear in time is however always present, and non-universal. 
The same exact term arises in the DP, and corresponds 
to the multiplicative contribution to the moments Z n = 

A 2 

e *™^-Rv(o) = e t™^fl 5 (o) arismg Ir0m t h e equa j rep ii ca 

(self-energy) terms after averaging the replicated parti- 
tion sum. Although usually dropped, these terms are 
present and correspond to an additive (non-universal) 

A 2 

non-random correction g-%i?^(0)t to \nZ. 
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